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Abstract. For non-equilibrium systems of interacting particles and for 
interacting diffusions in d dimensions, a novel fluctuation relation is derived. The 
theorem establishes a quantitative relation between the probabilities of observing 
two current values in different spatial directions. The result is a consequence 
of spatial symmetries of the microscopic dynamics, generalizing in this way the 
Gallavotti-Cohen fluctuation theorem related to the time-reversal symmetry. This 
new perspective opens up the possibility of direct experimental measurements of 
fluctuation relations of vectorial observables. 
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1. Introduction 

The understanding of non-equilibrium systems has been a main goal of statistical 
physics over the last decades, where a remarkable progress has been reached. Key to 
this advance has been the study of the statistics of several observables ranging from 
the entropy production or the work, to the density or the current (of particle, charge, 
energy, momentum...). As the probability distribution of these observables follow a 
large deviation principle, the study of their corresponding large deviation functions 
(LDFs) has become a subject of primary interest, because LDFs can be considered as 
the free energy analog in non-equilibrium systems [1,2]. For instance, the Macroscopic 
Fluctuation Theory [3], provides a variational formula from which one can compute 
the LDF of joint space-time densities and currents for driven diffusive systems, just by 
knowing two transport coefficients. 

Particularly important in the understanding of non-equilibrium large deviations are 
the so called fluctuations relations [4-15], valid arbitrarily far from equilibrium. Those 
relations basically exploit the time reversal symmetry of the microscopic dynamics, 
that gives rise to a relation between a positive fluctuation of certain observable and 
the corresponding negative fluctuation. However, for systems endowed with a spatial 
structure one may wonder if spatial symmetry yields relations between the probabilities 
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of vectorial observables in different directions. This was pointed out in [16], where an 
Isometric Fluctuation Relation (IFR) was obtained within the Macroscopic Fluctuation 
Theory framework under some assumptions and verified using numerical simulations. 
The IFR relation allows us to relate in a very simple manner any pair of isometric 
current fluctuations and it has been recently extended to the context of anisotropic 
systems [17]. In the IFR derivation the two dimensional system at hand was driven out 
of equilibrium by boundary reservoirs at different values in the presence of a held. In 
the present paper we consider systems driven far from equilibrium in the presence of a 
bulk held, i.e bulk-driven systems, and derive fluctuation theorems as a consequence of 
general spatial transformations, thereby extending the IFR case which was restricted 
to rotations. 

Fluctuations relations are very relevant because from them one can derive in the 
linear response regime the Onsager reciprocity relations and the Green-Kubo formulas 
and, even more important, other reciprocity relations beyond these two can be obtained 
by considering higher order response terms [19]. In a similar way the IFR relation 
implies a set of new hierarchies between current cunrulants and non-linear response 
coefficients [16]. 

The focus of this paper is on deriving spatial fluctuation theorems starting from 
the underlying stochastic microscopic dynamics. Whereas a clear understanding has 
been reached for the microscopic origin of the standard fluctuation theorem [5-8], it is 
an open problem to relate spatial fluctuation theorems to invariance properties of the 
microscopic dynamics. To achieve the result we will employ the Gibbs-formalism in 
space-time as introduced in [8]. The fluctuation symmetries which we derive here are 
a consequence of this space-time Gibbsian structure. They apply both in equilibrium 
context, when the space-time Hamiltonian is time-reversal invariant (see [20] for the 
analysis at equilibrium) as well as in systems driven away from equilibrium, for which 
time reversal symmetry is broken. 


2. Discrete systems 


We first consider particle systems embedded in 1 d following a Markov dynamics. For 
the sake of simplicity we restrict to discrete time, however similar results hold true for 
continuous time dynamics. The Markov process is denoted by {X(n) = (Xk(n))ken ■ 
n € N}, where X k (n) denotes the random position of the k th particle at time n. System 
configurations are denoted by the vector x = (# 1 , 2:2 , ■ • •) with the first particle located 
at site x\ € Kl = {1,..., L} d C if and similarly for all the other particles. During the 
evolution particles jump to their nearest neighbor sites with anisotropic probabilities. 
The anisotropy is tuned by the d-dimensional vector a = (a±,... ,ad) and in each 
direction we assume a weak asymmetry which is produced by a constant external field 
E = (Ei, ..., Ed). The dynamics is defined by the probability transition matrix with 


elements 


p(x,y) 


a s e ± £ if y = x k ' ±e “ 
0 otherwise, 


( 1 ) 


for k G N and s = 1,..., d. Here b(x , y) is a generic transition kernel from configuration 
x to configuration y , C is the normalization constant, e s is the unit vector in the 
s-direction and x k ’ dze ‘’ is the configuration that is obtained from the configuration x 
by moving the k th particle to Xk ± e s . 
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We consider an initial particle number (a conserved quantity of the dynamics) 
proportional to the volume, i.e. N = pL d , corresponding to a constant density 
0 < p < oo. We will be interested in the large deviations of the current vector 
Qm,l = z(N + ’ 1 — TV - ’ 1 , ..., N +,d — N~’ d ), where - ZV ±,S are the number of particle 
jumps in the ±s direction up to time M. This current is an additive functional 
of the Markov process (X(n)) ra6 N and it satisfies a large deviation principle that 
can be informally stated as P (Qm,l ~ Mq) ss e~ MIh ^ as M —> oo. By the 
Gartner-Ellis theorem [21] the large deviation function /^ can be obtained as 
Il(q) = sup A (A • q — pl( A)) where we define the scaled cumulant generating function 
of the current as 

Ml (A) := lim with Z m ,l{ X) =^(e x ' QM ’ L ), 


where E(-) denotes the average in path space. The application of Perron-Frobenius 
theorem allow us to express pl{ A) as the logarithm of the largest eigenvalue of the tilted 
matrix with elements p\(x,y ) = b(x,y)a s e ± (^ 7+ ^ if y = x k,±eB and p\{x,y) = 0 
otherwise. To study the system in the thermodynamic limit one needs to rescale by 
defining 


p{\) := 


lim 

L—foo 


Tl( A) 
L d ~ 2 


lim lim 

L—>o o M—>o o 


In Z m ,l( A) 
L d ~ 2 M 


( 2 ) 


Our first result provides the derivation of an anisotropic fluctuation relation 
in d dimensions from the assumption of a system invariance property. Consider N 
interacting particles whose dynamics is defined by the transition matrix (1). We denote 
by P the probability measure on the path space (i.e. the space of all trajectories) in 
the presence of an external field E = {E \,..., Ef), and Po the corresponding measure 
with E = 0. 


Theorem (Spatial fluctuation theorem). Let U : 7L d —> if he a transformation on 
the physical space. For a trajectory in path space x = (x(0), x(l),.. .x(M)), consider 
the bijective mapping : LI —> fl induced by U: {fl$x(n))k = Uxk(n ) . Assume that the 
current satisfies 

Qm,l(&x) = U Qm,l(x ) ( 3 ) 

and assume that Pq has the invariance property 


Po{x ) = P Q {&&) Vfgfi. 


( 4 ) 


Then the following fluctuation relations hold: VA 

Zm,l( A) = Z M , L (^{U~ 1 ) t (X + E)-E), (5) 

I L (q)-I L (Uq) = (U t E-E).q (6) 

where (t/ -1 )* denotes the transposed of the inverse. 

Notice that statement (5) is obtained at finite time and finite volume, whereas 
eq. (6) holds for any finite volume in the limit of large times. Furthermore, 
when the limit in (2) is finite, one has P(Qm,l ~ ML d ~ 2 q) « e~ MIjd 2/ ^ with 
I(q) = lim^^oo lL(q)/L d ~ 2 satisfying the analogous of eq. (6). We remark as well that 
to satisfy relation (4) the map U will also depend on the volume L and the anisotropy 
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a. However, to alleviate notation we do not write this dependence explicitly. In the 
proof we also shorthand Zm,l( A) as Z( A) and Qm,l(x) as Q(x). 

Proof. We start by observing that, from the previous definition of the current 
it holds that P(x) = exp [E ■ Q(x)]P 0 (x)/<yV, with jV = exp [E ■ Q(x)]P 0 (x). 

Applying this relation we have 

Z( A) = £ P{x)e X Q W = P °^l e (E+x) Q{s) . 


Thus, by using the invariance property (4) we get 


z(x) = £ 


Poi^x) 

jv 


0 (E+\)-Q(x) 


= £ P{@x)t 


-E-Q{Stx)+(E+X)-Q(x) 


Applying the change of variables y = Stfx and since & is a bijective map we find 
Z{X) = £ P(y)e _ - E ' Q( ^ +(B+A) ' Q (^” 1 ^. □ 
yen 


Hence, using the assumption (3), it is easy to check that (5) follows. By taking the 
limit M —> oo the same relation holds for y^ and therefore (6) follows by Legendre 
transform. 


3. Comments and examples 


If the transformation U is chosen as spatial inversion, i.e. Ui = — i for i £ Z d , then one 
recovers the standard Gallavotti-Cohen fluctuation relation, i.e. hl{X) = /j,l(—X — 2E). 
Notice that usually this relation is associated to time reversal invariance of the measure 
of the symmetric system. It was remarked in [7] that any transformation on path space 
such that M o — 1 would lead to the Gallavotti-Cohen fluctuation relation. The 
transformation on path space induced by spatial inversion in physical space has indeed 
such property. 

More generally, the theorem above allows us to deduce generalised fluctuation 
relations as a consequence of spatial symmetries, i.e. whenever a transformation U on 
the physical space satisfies (3) and (4) then (5) and (6) follow. To further illustrate 
this point we shall discuss examples of systems of non-interacting particles. In this 
case the dynamics can be studied in terms of a single particle and the scaled cumulant 
generating function can be explicitly solved, so that one can check by inspection 
which spatial symmetries hold. For instance, for a system of independent random 
walkers (RW) where each particle at site i jumps to site i ± e s with probability a s e ± ~E 
with periodic boundary conditions, an elementary application of Perron-Frobenius 
theorem gives 


y™(\) = pL‘ i In 


a 

n 

,s=l 


a,e 


(E b +\ s )/L 


+ a s e 


~(E S +X S )/L 


For d = 2 and by doing a change of variables to polar coordinates (z, 9) such that 
Xi = zcosO — Ei and X 2 = zsin9^/a — E 2 , with a = ai/a 2 being the anisotropy ratio, 
the above expression reads 


/i^ w (z, 9) = pL 2 In «i ‘ + e L j 


+a 2 




(7) 



Spatial fluctuation theorem 


5 


In the isotropic case (a = 1) one recognizes by inspection the discrete symmetries 
of the system leading to a fluctuation theorem. Namely, /i£ W (z,0) = 0'), 

for 0' = m7r/2 zb 6 , Vra G Z (see Fig. 1). It is then natural to expect that 




n RW / 

Figure 1 . (Color online) Plot of d _ 2 


in d = 2 as a function of 0 = 


arctan(-^= )• In both figures p = 0.5, E = (5,0), z = 4 and the anisotropy 

ratio is a = 1 in the top and a = 2 in the bottom. The inset of the top panel 
shows the 8 discrete symmetries of the system in A-space for L = 10, starting 
from an initial angle e.g. 6 = 7r/6, corresponding to the transformations U 
given by the diagonal and the anti diagonal matrices with elements ±1. As L 
increases the discrete symmetry scales to a continuous symmetry associated to the 
cos 6 — yfa sin 6 

4= sin 0 cos 0 


transformation U = 


and (9) is then satisfied. 


for diffusions one is lead to consider continuous spatial symmetries. This can be 
seen by considering the diffusive scaling limit of the previous example. If we define 
{X( L \n) = (Xl L \n))ken ’■ n £ N} the positions at time n of the N independent 
random walkers (labeled by k ) in a volume of linear size L, then the process 

{R t = (Rk,t)ke n : t G R + } defined by R t := lim^oo A ' ( ^ L will be a family 
of N independent anisotropic Brownian motions satisfying the stochastic differential 
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equation dRk.^t — AEdt + \/~AdBk t t with A the d x d diagonal diffusion matrix with 
elements A s>s = 2 a s and B^ t £ denoting independent standard Brownian motions. 
An immediate computation gives 

lim dL d _ 2 = H BM (X) = pVa s A s (A s + 2 E s ) (8) 

L—foo L a z z ' 

S — 1 

where /i BM (A) := limr^,*, In E(e A '^ T ) , being Q T = J2k=i Rk,T the current up 
to time T. From the explicit expression (8) we see that a spatial fluctuation relation 
holds, i.e. 

H BM {\) = f i BM ((U- 1 ) t (\ + E)-E^ (9) 

with U such that UAU t = A. For d = 2, U takes the form given at the end of the caption 
of Fig. 1. As we shall see below, such a relation can be traced back to the invariance of 
the path measure of the anisotropic Brownian motion under a spatial transformation 
U such that UAU 4 = A ,in agreement with the findings of [16,17] in the context of 
the Macroscopic Fluctuation Theory [3]. Notice indeed that the MFT predicts for 
a two dimensional anisotropic periodic system in the presence of a hied E that the 
probability of having a time-integrated current fluctuation qo = ^ J 0 T dt f drq(r , t ) up 
to a long time T is P{qo) ~ e -TL I ( q A with the LDF [17] (see also [18] for the one 
dimensional case) 

J(g 0 ) = lim inf ^(p,g), (10) 

T—>-oo (p,g)GO 1 

where 

^(p, < i) = \j dt j dr ( q + d (p)Vp ~ - Ecr (d)) T(J (p) _1 (9 + d (p)^p - E°{p))- 

Here the diffusivity and the mobility are given by D{p) and cr(p) respectively. In (10), 
is the set of paths (p, q) satisfying the continuity equation and whose time-integrated 
current is qo, 

V = | (P, Q) ■ ^ dt 

The cumulant generating function is then given by the Legendre transform of (10) 

p(X) = sup[Ag 0 - I(q 0 )]. (11) 

Qo 

In particular, for N anisotropic independent random walkers in a volume L d such 
that p = N/L d , D{p) and a(p) are given by the diagonal matrices with elements a s 
and 2 pa s respectively. In this case, it is easy to check that by solving the variational 
problem given by (10) we get (8) as expected. 

4. Diffusions 

To state a (spatial) fluctuation relation for systems following a generic diffusion process 
we consider an abstract path space ft, and a bijective measurable transformation 


J drq(r, t) = q 0 , d t p = -V • g j 
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Si : Q —> with inverse Si~ x . Elements in the path space are denoted by oj. For 
diffusions we will consider fl = ^([0, T], R d ), i.e. the set of all continuous paths up 
to time T taking value in IR d , and (S$(u>))t = Uuit, with U : an invertible 

spatial transformation. 

Proposition. Consider a probability measure on path space Q of the form P(dui) = 
e H (y>)P 0 (du>) where Po is Si-invariant. For all ip : fl —>■ K. we have the identity 

E(e v °^) = E(e v e ffo ^ 1_ff ) ( 12 ) 

Proof. 

J e v> ° m dP = J e tpo ^e H dP 0 = J e^e^^dPo 
= J e v> e Hom ~ 1 ~ H dP. □ 

In the particular case that tp = —7(17 o Si~ x — H) with 7 G I, the relation (12) 
gives E(e 7 ^°^^ ff ^) = E(e^ 1_7 ^' tfo ^ ~ H ^). Even more, specifying to a transformation 

such that Si = Si~ x this identity is exactly a symmetry of the form of the standard 
fluctuation theorem for the quantity H o Si — H. 

Denoting by Po J 1 the image measure of P under Si, i.e., J fd(P o&) = f (f o S#)dP, 
it can be readily verified that dP = e H dP 0 implies that d(P o Si) = e Hom 1 dP 0 . 
Therefore, as we shall use below, we can write 

Hogft- 1 ~h _ d(P 0 Si) _ d(P o Si) / dP 

dP ~ dP 0 / dP 0 ' [ ’ 

The abstract setting of the proposition above can be used to derive the spatial 
fluctuation theorem for finite time T in the context of interacting diffusions. We shall 
illustrate this by considering the overdamped Langevin dynamics {X t = (X k , t ) ke N : 
t € R + } describing N particles (labeled by k ) subject to a drift vector which can arise 
from an applied force to each particle F and/or a conservative potential V, with a 
positive definite constant diffusion matrix A. Then the stochastic differential equation 
for the k th particle reads (Ito convention) 

dX k>t = F(X kJ )dt + V k V(X t )dt + VMB ktt (14) 

where B kt is again a standard Brownian motion. Notice that V can model a self¬ 
potential as well as an interaction potential. To obtain the new process applying a 
transformation Si of the type (Si(ui))t = Uuit, we remind that if X is multivariate 
normal distributed with mean zero and covariance matrix A, then for any d x d matrix 
U, UX is multivariate normal distributed with mean zero and covariance UAU t . As a 
consequence if Y k t = IJX k t then 

dY k<t = UF(U~ 1 Y k j)dt + US7 k V (U~ 1 Y t )dt + VuMPdB ktt 

where we denote U~ 1 Y t the collection U~ 1 Y k ^ t . Vfc. For the process Y t to be absolutely 
continuous w.r.t. the X t process, we need that UAU* = A. Moreover, assuming that 
the potential is invariant under the transformation 17, i.e. V(Ux) = V(x), the process 
Yfc.t satisfies 


dY k , t = UF(U~ 1 Y ktt )dt + V k V(Y t )dt + y/AdB M . 


( 15 ) 
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Thus, the process whose paths’ distribution is invariant under the transformation £% is 

dZkj = S7kV(Z t )dt + VAdBk,t (16) 

By using the Girsanov formula [22] one can compute dPx = e H dPz , i.e. the relative 
density between the path space measure Px of the process (14) and the path space 
measure Pz of the process (16). Analogously the measure Py of the process (15) and 
the measure Pz are related by d(P x ot%) = dPy = e HoSi 1 dPz- Hence, if we denote U 
and A as the N x (d x d) block diagonal matrices with elements U and A respectively 
on the diagonal and F as the N x d column vector consisting in copying N times the 
vector F, then by applying (13) we get 

(HoAg~ l -H)(u) = [ A~ 1 (UF(U~ 1 cj t ) — F(co t ))duj t 

Jo 

A-\UF{U-^ t ) - F(u t )) ■ W (u t )dt 

T 

(i UF{U~ l uj t ) • A~ 1 (UF(tj~ 1 u t )dt 

+ \ f F(u t ) ■ A~ 1 F(uj t )dt. (17) 

* Jo 

In addition, if the force F is constant we find (we put wo = 0) 

(H o &- 1 - H)(u) = A~\UF - F) ■ Q t (oj) (18) 

with Qt(oj) = )CfeLi — / q T VkV(<jJt)dt^J. Furthermore, by choosing tp(ui) = 

A • Qt(ui), with A e R d , we have that (ip o &)(uj) = A • UQt(oj) = U t X ■ Qt{oj) 
and from (12) and (18) we get E(e c/ A < ? T ) = E(e^ A+A (uf—f)-q t ^ gy defining 
Zt{ A) := E(e A '^ T ) it readily follows that 

Z T ( A) = Zr^U-'Yi A + A~ 1 F) - A~ 1 F^j (19) 

which is the analogous relation to the previously found for the discrete setting (5). 
Notice that if A is the diagonal matrix with elements A s s = 2 a s , F = AE, V = 0 and 
/i BM (A) := lim 7 ’_ >00 (TL <i ) _1 lnZ^A) we recover from the above equation the previous 
result (9) for the anisotropic Brownian motion. 

5. Conclusions 

In this work we have derived a spatial fluctuation theorem (SFT) for interacting particle 
systems and interacting diffusions driven out of equilibrium by an external field. It has 
been proved in both cases that the SFT can be traced back to an invariance property 
of the microscopic path space measure under spatial transformations. Remarkably, this 
result holds for finite time in the case of interacting diffusions (19), with UAU 4 = A 
and V(Ux) = V(x), as well as for finite volume for interacting particle systems (5). 
In the latter, the spatial transformations yielding the SFT correspond to discrete 
symmetries associated with the underlying lattice geometry of the system at hand. As 
the system linear size increases these transformations become continuous. 

The SFT gives a new perspective on how the microscopic symmetries of the system 
are reflected at the fluctuating level. Whereas the standard fluctuation theorem is 
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based on time reversal-symmetry of the microscopic dynamics, the SFT proves that 
the spatial symmetries have also a word to say. It is worth emphasizing that from the 
SFT new hierarchies for the currents cumulants and for the non-linear response can be 
obtained [16,17]. For more details see section 6.2 of the J. Stat. Phys. paper of [16], 
where these hierarchies are derived from the SFT, /z(A) = n^(U~ 1 ) t (\ +E) — u) , when 
the transformation U is a rotation. To obtain the hierarchical cumulants relations 
that has to consider the limit of infinitesimal rotations, namely U = I + A6LSf, with 
I the identity matrix and Jz? any generator of d-dimensional rotations. For instance, 
to the lowest order these hierarchies imply Onsager’s reciprocity symmetries and 
Green-Kubo relations for the linear response coefficients, with the additional prediction 
that, the linear response matrix is in fact proportional to the identity. In addition, the 
bare extension to more than one dimension of the standard fluctuation theorem only 
allows to compare one spatial direction and its reverse. It is within the SFT that one 
can relate probabilities of observing two current values in arbitrary different spatial 
directions (see [23] for the experimental measurements of the current fluctuations of 
a self-propelled rod immersed in a sea of spherical beads). From the experimental 
point of view, as the standard fluctuation relation leads to the the study of free-energy 
differences in terms of work distributions (see e.g. [24,25]), one can expect that SFT 
contains further information on the distribution of the mechanical work. The study of 
the full Langevin dynamics (including inertia) still remains as an open problem, as 
well as systems subject to noise that is non-Gaussian. 
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